Current quark mass effects on chiral phase transition 
of QCD in the improved ladder approximation 
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Current quark mass effects on the chiral phase transition of QCD is studied in the improved 
ladder approximation. An infrared behavior of the gluon propagator is modified in terms of an 
effective running coupling. The analysis is based on a composite operator formalism and a varia- 
tional approach. We use the Schwinger-Dyson equation to give a "normalization condition" for the 
Cornwall-Jackiw-Tomboulis effective potential and to isolate the ultraviolet divergence which ap- 
pears in an expression for the quark-antiquark condensate. We study the current quark mass effects 
on the order parameter at zero temperature and density. We then calculate the effective potential 
at finite temperature and density and investigate the current quark mass effects on the chiral phase 
transition. We find a smooth crossover for T > 0, fi = and a first-order phase transition for fi > 0, 
T — 0. Critical exponents are also studied and our model gives the classical mean-field values. We 
also study the temperature dependence of masses of scalar and pseudoscalar bosons. A critical end 
point in the T-fx plane is found at T ~ fOO MeV, y. ~ 300 MeV. 
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I. INTRODUCTION 



Recently there has been great interest in studying the phase structure of quantum chrome-dynamics (QCD). We 
expect that at sufficiently high temperature and/or density the QCD vacuum changes into a chirally symmetric/ 
deconfinement phase jl| , a color superconducting phase (^-^] or a color-flavor locked phase || . They may be realized 
in high-energy heavy- ion collisions at the BNL Relativistic Heavy Ion Collider and CERN Large Hadron Collider. 
These phase transitions are important also in the physics of neutron stars and the early universe. 

The phase structure of QCD has been studied by various methods. At finite temperature the lattice simulation is 
Qh! powerful for studying the phase structure of QCD. On the other hand, the so-called QCD-like theories, one category 
^ of the effective theories of QCD, are still useful for studying the chiral phase structure at high temperature and/or 
density [6-15]. QCD in the weak coupling limit is utilized to study the color superconductivity [2-5]. However, most 
studies have been done in the zero external field limit; i.e., in the chiral limit for the case of the dynamical chiral 
symmetry breaking. It is, then, desirable to study a more realistic situation where the chiral symmetry is explicitly 
broken by a current quark mass. In this paper we neglect the quark pairing and study the current quark mass effects 
on the chiral phase transition between SU(N C ) x SU(Nf)i, x SU(Nf)n and SU(N C ) x SU(Nf)L + R. In order to study 
nonperturbative phenomena such as the dynamical chiral symmetry breaking, the Schwinger-Dyson equation (SDE) 
or the effective potential for a composite operator has been widely used. However, in order to find the true vacuum 
it is necessary to calculate the effective potential. Furthermore, in the studies of the SDE, it was known that there 
is a difficulty in removing a perturbative contribution which is quadratically divergent from the order parameter of 
chiral symmetry Jl6| , p7[ . 

Our analysis starts from the Cornwall-Jackiw-Tomboulis (CJT) effective action for a composite operator Jul and 
QCD in the improved ladder (rainbow) approximation [19-21]. The improved ladder approximation is the simplest 
nonperturbative truncation scheme which is consistent with the renormalization group. At zero temperature and 
density, it reproduces the physical quantities which are insensitive to the model parameter p2| . Therefore, we expect 
that it is valid also at finite temperature and/or density. We use the SDE to modify the CJT effective potential. 
The result is used to give a "normalization condition" for the effective potential [E3] and to isolate the ultraviolet 
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divergence which appears in an expression for quark-antiquark condensate. The divergence depends solely on the 
qurak mass and is independent of the condensate. 

This paper is organized as follows. In Sec. II we derive the modified form of the CJT effective potential for quark 
propagator at zero temperature and density. The critical value of the coupling constant and the critical number of 
quark flavors are examined for comparison with the results by other authors. We then extend the modified effective 
potential to finite temperature and density. In Sec. Ill we first determine the value of A QCD by a condition f„ = 93 
MeV at T = /x = and in the chiral limit. Then we calculate the effective potential numerically and investigate 
the phase structure for nonzero T, fi and nonzero current quark mass. We also examine the critical exponents, the 
temperature dependence of masses of the scalar and pseudoscalar bosons and the critical end point. Section IV is 
devoted to the summary and discussion. We fix the mass scale by the condition A QCD = 1 except for numerical 
calculations. 



II. EFFECTIVE POTENTIAL FOR THE QUARK PROPAGATOR 

A. CJT effective potential at zero temperature and density 

At zero temperature and zero density, the CJT effective potential for QCD in the improved ladder approximation 
is expressed as a functional of S (p) the quark full propagator jl9| : 

V[S] =V 1 [S] + V 2 [S], (1) 

Vx[S] = I ^ Tr{\n[S Q 1 (p)S(p)} - S Q 1 (p)S(p) + l} , (2) 

V 2 [S]=- l -C 2 J J 1 0^ 1 0^f( p , q )Tr[ 1 ,S{p) lv S{ q )]D^{p~ q ), (3) 

where C 2 — (N 2 — 1)/(2N C ) is the quadratic Casimir operator for the color SU(N C ) group, So(p) is the bare quark 
propagator, g 2 (p, q) is the QCD running coupling of one-loop order, D fJ,u (p) is the gluon propagator which is diagonal 
in the color space and, "Tr" refers to Dirac, flavor, and color matrices. The two-loop potential V 2 is given by the 
vacuum graph of the fermion one-loop diagram with one-gluon exchange. 
After Wick rotation, we use the following approximation 

g 2 (PE,qE) = 0(p E - qE)g 2 (pE) + 0{q E -pE)g 2 (qE)- (4) 



This approximation is occasionally called Higashijima-Miransky approximation |19|J20[| . In this approximation and 
in the Landau gauge, renormalization of the quark wave function is unnecessary and the CJT effective potential is 
expressed in terms of T,(p E ), the dynamical mass function of a quark as follows 



V[£{p E )} = V x [E(p B )] + V 2 [£(p E )], (5) 
^i[S(p B )] = -2 J 

1 W 4 -n„ V.(<n„MY,t<n„\ - <m( W 

(6) 



V 2 [£(p E )] = -6C 2 J J 



d 4 PE . ^ 2 (pe)+Pe 
(2tt) 4 n m 2 (A)+p| 

A ^pe S(p £ )[S(p g )-m(A)] 

(2^)4 ^(pe)+P 2 e ' 

A ' A d*p E d*q E g 2 ( PE ,q E ) 



(27r)4(27r)4 {p E -q E f 



^( Pe )+P 2 e m<lE) + q 2 E 



(7) 



where an overall factor (the number of light quarks times the number of colors) is omitted and m(A) is the bare quark 
mass. In the above equations we temporarily introduced the ultraviolet cutoff A in order to make the bare quark 
mass well defined. 

The extremum condition for V with respect to the variation of Y,{jp E ) leads to the following SDE for the quark 
self-energy 

w \ (\\ , f A di lE 9 2 {PE,qE) Z(q E ) 

Z(PE) =TO(A) + 3C 2 / vT 7 yT \ j 2~ ' 8 

J 2tt) 4 {p E - q E ) 2 £ 2 + q E 
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In the Higashijima-Miransky approximation, since the argument of the running coupling has no angle dependence, 
we first perform the angle integration using the result: 



^4 n , s 1 nl \ 1 

7 = 0(p E - qE)— + 0{q E - Pe) — , 

[PE - qE) 2 Pe Te 



where 



2tt 2 



dip sin ip / d9 sin 9 I dp. 



Then Eq. (||) is reduced to the following differential equation : 

E(p B ) _ (4tt) 2 d / 1 d£(p B ) 



£ 2 fe)+p| 3C 2 p\d v \ \A(p E ) dp% 
with the two boundary conditions 

1 dH{ PE ) 



&(pe) dp\ 

V{p E ) dX( PE ) 
Mpe) dp 2 E 



0. 



PE=0 



= m(A), 



PE=A 



where the functions 



and 



V{ PE ) 



A(p B ) 



9 2 (pe) 
Pe 



dp 



2 v (pe), 



are introduced. 

Substituting Eqs. (|) and |[l]) into Eqs. (|) and (@), we obtain 



V[E(pi 



-2 
+2 



A dV , s 2 ( P£ )+p| 

(2tt) 4 n m 2 (A)+p 2 E 
A d*p B E(p B )[E(p B )-m(A)] 



(27T) 



Z 2 (pe)+P 2 e 



A dV ln s 2 ( P£ )+p| 



(2tt) 4 m 2 (A)+p| 



2(4tt) 



3C 2 



— ^ - m(A)] ^4^2" 



(2^) 4 



p|dp| V A (Pis) d P : 



= -2 



A d 4 PE Z 2 ( PE )+p 2 E 
(2tt) 4 n m 2 (A)+p|; 



2 



A 



1 



where we used a partial integration in the last line and 



1 d^( PE ) 

2 
E 



V s = F(A) - F(0), 

F( Pg ) = ^ r [S( Pg )-m(A)] 1 y . 
3(72 dp E 



(9) 
(10) 

(11) 

(12) 
(13) 

(14) 
(15) 



(16) 

(17) 
(18) 
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Hereafter, we consider the effective potential in the continuum limit (A 
using the running coupling 



f(PE) = 



2n 2 a 



24 



and the corresponding asymptotic form of the mass function 



2n f 



oo). Let us begin by evaluating F(A) 



(19) 



S(p B ) -► m(A) 



/lnp 



VlnA 2 



-a/2 



a ft 2 

Pe 



2 \a/2-l 



(20) 



where rif is the number of quark flavors which contributes to the running coupling. Throughout this paper, with the 
exception of discussion on a critical number of massless flavors, we put N c = rif = 3, namely a = 8/9. As we will 
show below, the parameter a is related to the order parameter of chiral symmetry a.k.a. the quark condensate. Note 
that Eq. (|2(]) is to be understood in the sense that for exact chiral symmetry the first term on the RHS is zero and 
the second term is the dominant one, while in the presence of explicit chiral symmetry breaking the first term gives 
the dominant asymptotic behavior. With the explicit chiral symmetry breaking, there will be many terms which are 



suppressed by 
Using Eqs. ( 



powers of ln(p E ) compared to the first term, but dominate the second term as p E — > oo Jit 



19D and (|2C|), we obtain 

2 
3Ch 



F(A) 



A 2 1 ' lnA 2 + l 



27r 2 a 



a 
2 



>(A) 



A2hiVJ -^(lnAY /2 - 2 (lnA 2 + l-a/2) 



67r 2 aC 2 



A 2 (lnA 2 ) a / 2 



a m(A) ^ (lnA 2W2-i 
2 A 2 In A 2 A 4 ' 



(21) 



We note that m(A) is the bare quark mass defined at the scale A, as mentioned before, and the factor (In A 2 ) a / 2 m(A) 
in Eq. ( pl| ) is equal to (In K 2 ) a / 2 mp ! {K) which is cutoff independent. [] Hence, F(A) vanishes in the continuum limit, 
provided a remains finite in this limit. As concerns F(0), it turns out to be 



F(0) 



3C. 



\E(p E ) - m(A)] 
(Pi) 2 



dH(p E ) 



P E d 9 2 (PE)/dp 2 E - g 2 {p L 



d P 2 E 



PE=0 



(22) 

and (|35|), we can set 



Since we will introduce an infrared finite running coupling and mass function in Eqs. ( 
F(0) = 0. 

After all, in the continuum limit, we get Vs = and the modified version of the CJT effective potential becomes 



V[E(p f 



d A PE Z 2 ( PE )+p 2 E 
(2tt)4 



Pe 



2 f 1 



dp E 



(23) 



1 In this paper, we use the mass-independent renormalization scheme |^,^5). In this scheme, all the renormalization constants 
are fixed by massless theory. The bare quark mass m(A) and the renormalized mass ?ti_r(k) at the renormalization point k are 
related as 

m(A) = Zg 1 (A, K)mn(K) 
where Zg 1 (A, k) is the renormalization constant for the composite operator qq: 

(qq)n = Zg 1 ■ (qq)A, 

Zs\A,k) 



In k 
In A 2 



2 \ a / 2 
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A few comments are in order. 

(1) The extremum condition for U[£(p£;)] in Eq. with respect to T,(pe) leads to Eq. (|li"| ) which is equivalent 
to the original equation (||) in the Higashijima-Miransky approximation apart from the two boundary conditions. We 
will take account of these conditions when we introduce the trial mass function. 

(2) Even if chiral symmetryis explicitly broken, we do not use the condition adopted in Ref. p3[ . Instead, our 
condition is V s = (see Eq. (fL7|)). 

The relation between a and the quark condensate is as follows. The renormalization group invariant vacuum 
expectation value of qq is given by 

(qq) = (0\(qq)R\0)(\nK 2 r^ 2 

= lim (0|(^) A |0)(lnA 2 )- a / 2 

A— »oo 

= lim G(A), (24) 

A^oc 



where G(A) is defined as 



with 

Z(A) = (lnA 2 ) Q / 2 . (26) 

Using Eq. (|ll|), G(A) is rewritten as 

G(A, . - Z -(A)lf f A ' *& < (>jSfO) . (27) 
3G 2 Jo d ?4 V a (pb) ^ / 

This expression is convenient to isolate the ultraviolet divergence as the RHS is linear in £(ps)- Using Eqs. ( |l^ ) and 
(pp|), as can be verified by direct calculation, G(A) is obtained as 

G(A) = m(A)/ p (A) + cI„(A), (28) 

where the lower bound in the definite integral automatically vanishes by the same argument as given for F(0) and a 
perturbative contribution I p (A) and a nonperturbative one J n (A) are given as follows: 



^(A) = -^^A 2 (lnA 2 )- a / 2 , (29) 

IN 
3ir z aC2 

As discussed in Ref. , the perturbative contribution is quadratically divergent. We are interested in the purely 

nonpertubetive effect; therefore, we simply subtract the perturbative effect from Eq. (|2|) and re-define the quark 
condensate as 



(qq) = aI n (A) 
3 



2ir 2 a 



(31) 



The relation between a and (qq) is the same as in the massless theory. We do not use Eq. (^) but Eq. (|3l]) to 
extract (qq) from a m i n , the location of a minimum of the effective potential. Consequently, the asymptotic behavior 
of the mass function in the continuum limit takes the following form 



nPE) - mR (\npl)- a / 2 - (qq)^(lnplr^\ (32) 



where 



m R = m R (n){\nK 2 ) a ' 2 = lim m(A)(ln A 2 )"/ 2 (33) 

A^oo 
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is the renormalization group invariant quark mass. This asymptotic behavior is completely coincident with that 
obtained by using the operator product expansion and the renormalization group equation [[26| . Notice that in these 
arguments there are no ambiguities in separating the perturbative and the nonperturbative contribution in (qq). 

Now we are in a position to introduce a modified running coupling and a trial mass function. We use the following 
modified running coupling El]] 



2n 2 a 



f(PE) = TTT^T^v (34) 



where p R is a parameter to regularize the divergence of the QCD running coupling at pe = 1(A QCD ). This modified 
running coupling approximately develops according to the QCD renormalization group equation of one-loop order at 
large p E , while it smoothly approaches a constant as p 2 E decreases. 

Corresponding to the above running coupling, the SDE with the two boundary conditions suggests the following 
trial mass function 

T,(p E ) = mnlHpl+pl)}-"/ 2 + ^-^l H p E +pl)}^-\ (35) 

Pe "t Pr 



The trial mass function is infrared finite, and moreover, it has the same asymptotics as in Eq. (|32^ 

When we also consider the pseudoscalar degree of freedom in the effective potential, we expand the full quark 
propagator in Minkowski space as follows 

iS^ip) = 7 -P- ^(p) - ij 5 E 5 (p). (36) 

An effective potential for this case is obtained straightforwardly as 

d A p E , ^ 2 (pe) + ^ 5 (pe)+P 2 e 



V[£(p E ),£ 5 (p £ )] = -2 Jj^ In- 



Pi 



-3C 2 J dpE Afe) 



(37) 



The mass function Y,^(j>e) satisfies the SDE which has the form as in Eq. ( |ll| ) where the bare quark mass is set zero. 
Therefore the trial mass function for Y,^(p E ) should be 

Z 5 (pe) = -^^[Hpe+pI)T' 2 -\ (38) 
Pe ~t~Pr 

where 175 = 2ir 2 a(qij 5 q) /3. However, except when we consider a mass of pseudoscalar boson we can put T, 5 (pe) = 
without loss of generality. 

Substituting Eqs. (|34|) and (|3^) into Eq. (p3|), we have the following expression for the effective potential 
V(a;m R ) = —7—^ dt e*(e* - e tR ) In ■ 



(47r) 2 J tjt e* - e tR + m 2 R t- a 

m R a f°° e-*(e* - e*«) 2 (f + 1 - a/2) 



At: 2 J tR ' t{e l - e*« + te l 

9 2 r ^ (l-e'«-') 2 (t + l- a /2) 



/ dt 5 — , / , rr^-, (39) 

32?r 2 J tR t 2 - a (e l -e'« +te t ) 

where we made the change of the integration variable 

t = \n(p 2 E +p 2 R ) , t R = \n(p 2 R ), (40) 

and subtracted the a independent term which is exponentially divergent in t. The infrared regularization parameter 
t R specifies the coupling constant in the low energy region. Indeed, g 2 (Q) = 2Tr 2 a/t R is proportional to the inverse 
of tR. We note that in the first two terms in Eq. (^) there are logarithmic divergences in t which cancel out each 
other. This is because the fact that we take into account the correct renormalization group effects in s 2 (pe) and in 

Before we turn to the effective potential at finite temperature and density we examine some numerical consequences 
of Eq. (||^) and make a comparison with the results by other authors. Figure 1 shows the t R dependence of (qq) in 



G 



the chiral limit. Chiral symmetry is broken for relatively small tR whereas it is restored for large tR. Notice that (qq) 
is stable under the change of tR if tR < 0.4 and the critical value of tR is {tn)critical — 0.56. This critical value of tn 
is about one third of that obtained in Ref. |2l| in which the SDE with the same running coupling has been solved. In 
our model, 

2n 2 a 

\&s jcritical 



critical 

= 3 w < 41 > 

where 7r/(3C2) is the critical value obtained from the SDE analysis in Ref. J2l| and also in quenched QED 4 with the 
replacement a by C^a p6|. We note that, in Ref. |19 2lJ], it has been shown that the {a s ) cr i t i C ai obtained by using 
the CJT effective potential (||) and (Q) is about two times of w/ (3C2). These differences may arise from the fact that, 
in the framework of the variational approach, we restricted the class of trial function and we used the modified form 
of the CJT effective potential. 

Figure 2 shows the plot of the quark condensate (qq) versus quark mass wir analogous to the relation between 
the spontaneous magnetization of the ferro-magnetism and the external magnetic field. The curves show the cases 
tR = 0.1 and tR = 0.8. When tR is sufficiently small, i.e., the coupling constant in low energy region is sufficiently 
large, the quark condensate discontinuously changes its sign at wir = 0. This is a clear evidence of the dynamical 
chiral symmetry breaking. On the other hand, when tR is sufficiently large, the discontinuity of the quark condensate 
disappears. 

The chiral phase, even at zero temperature and density, depends on the number of massless flavors. At some value 
of nj (less than HA c /2), a phase transition to the chirally symmetric phase is expected. In Ref. |27j], it has been 
argued that, in connection with the infrared fixed point in the two-loop f3 function, there exists a critical number of 
fermions n c j , above which there is no chiral symmetry breaking. Figure 3 shows the tR dependence of n °j , where n ^ is 
considered as a continuous number. For tR = 0.1 our model suggests that nfj = 10, and this result is consistent with 
lattice QCD results p§ which gives 8 < n c f < 12 (see also Ref. p7|]). 



B. Effective potential at finite temperature and density 

In this subsection we discuss the effective potential at finite temperature and density. In order to calculate the 
effective potential at finite temperature and density we apply the imaginary time formalism |29| 



n— — o 



/K + i/i), (42) 



where u> n = (2n + l)irT (n = 0, ±1, ±2, • ■ •) is the fermion Matsubara frequency and /i represents the quark chemical 
potential. In addition, we need to define the running coupling and the trial mass function at finite temperature and 
density. We adopt the following functions for T>T,fj,(p), Tit^(p) and Ti^t^(p) by replacing P4 in V(pe), ~E(pe) and 
^h{pe) with Lu n : 



V tAp) = ^,,,2 , 1^2 , -2^ ..2 , M 2 ' ( 43 ) 



2n 2 a 1 

^ T Ap)^^RH^+\p\ 2 +p 2 R)v a/2 



[Huj 2 n + \p\ 2 +pl)] a/2 -\ (44) 

n, 

\P\ 2 + P 2 R 



u>l + \p\ 2 +p 2 R 

^tAp) - ttSt ~ M«l + W +P 2 r)T' 2 ~ 1 ■ (45) 



In Eq. ( f43| ) we do not introduce the /i dependence in T>tAp)- The gluon momentum squared is the most natural 
argument of the running coupling at zero temperature and density, in the light of the chiral Ward-Takahashi identity 
[^,|lj. Then it is reasonable to assume that T>tAp) does not depend on the quark chemical potential. In addition, 
the screening mass is not included in Eq. (|43|). We comment on this point in Sec. IV. 

As concerns the mass function, we use the same function as Eqs. ( |35| ) and (^) except that we replace p^ with w„. 
As already noted in Sec. II A, the quark wave function does not suffer the renormalization in the Landau gauge for 
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T = fi = 0, while, the same does not hold for finite temperature and/or density. However, we assume that the wave 
function renormalization is not required even at finite temperature and/or density, for simplicity. 

Furthermore, we neglect the T-\i dependent terms in the quark and gluon propagators which arise from the per- 
turbative expansion. We expect that the phase structure is not so affected by these approximations. 

Using Eqs. (]43|), (Q) and it is easy to write down the effective potential V(a, 05; ran) (see Appendix). Even 
at finite temperature and/or density, we determine the value of (qq) through (qq) = — (3/2ir 2 a)a m i n where a m i n is 
the location of the minimum of V(a, 05; ttir). 



III. CHIRAL PHASE TRANSITION AT HIGH TEMPERATURE AND DENSITY: NUMERICAL 

RESULTS 



In numerical calculation, as mentioned before, we put N c — rif = 3. Furthermore, since it was known that the 
temperature and chemical potential dependence of the quantities such as (qq) and f n are stable under the change of 
the infrared regularization parameter If |, in the first place, we fix t R = 0.1 and determine the value of A QCD by the 
condition f v — 93 MeV at T = fi = and m R — 0. In this case, the pion decay constant is approximately given by 

d 4 p E S(pb) / w x Pe dT,{p E ) 



and we have A QCD = 738 MeV. Secondly, we assume the light quarks (u and d) are degenerate in mass and take the 
quark mass evaluated at n = 1 GeV as 

, m u (lGeV) +TO d (IGeV) 
m fl (IGeV) = — '— = 7MeV. (47) 

With Eq. ( |33"|) the renormalization group invariant quark mass ma extracted from the above-mentioned value becomes 

m R = 7.6 x I(T 3 A QCD . (48) 



A. T / 0, n = case 



The phase diagram in the chiral limit is shown in Fig. 4. We have found second-order phase transition, with the 
critical temperature T c = 129 MeV for fi = Jf|]. 

Figure 5 shows the temperature dependence of the effective potential at \i = for mij(IGeV) = 7 MeV. We 
can realize that a m i n the minimum of the effective potential continuously approaches to the origin as temperature 
increases. Figure 6 shows the temperature dependence of a m i n for the chiral limit and the case mu(\ GeV) = 7 MeV. 
Below T c , the two curves are almost the same apart from the small difference in their magnitude. On the other hand, 
above T c , while er m i n remains nonzero for mfl(IGeV) = 7 MeV, the one in the chiral limit vanishes. In other words, the 
temperature dependence of the order parameter is dominated by the dynamical symmetry breaking in low-T region, 
whereas by explicit breaking in high-T region. 

Let us examine the critical exponents. Assuming the mean-field expansion, we expand V as follows 

V = a 2 (T)<r 2 + a 4 (T)a 4 + b(T)m R a + ■■■, (49) 

where ellipsis represents higher order contributions in m R and a that is expected to be small near T c and the critical 
temperature is determined by the condition a 2 (T c ) = for 777,7? = 0. The three critical exponents (3, 7 and 5 are 
defined as follows 



T 







(qq)\m R =o^ , (50) 



9{qq) 



drriR 



(qq)\ T=Tc ~m R /5 , (52) 
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where T < T c . Since we have confirmed that, except for ay (T), all the coefficients in Eq. (^) do not have singular 
behavior near T c we can define the exponent 7, instead of Eq. (|5lJ), as 



T 



7 



a 2 (T)^ (^1 — J . (53) 

These critical exponents are determined numerically by the use of the linear log fit, for instance, 

ln(qq) | roR=0 = pin(l-^-)+C, (54) 

where C is independent of T. In order to determine the exponent we use the \ 2 fitting (see Figs. 7 and 8). The 
results are 

p ~ 0.5 , 7 ~ 1 , (5-3. (55) 

Therefore, our results confirm the Landau theory of the second-order phase transition, i.e., a second-order phase 
transition with the classical mean-field values for the critical exponents, and coincident with that obtained in Ref. 
II- 

In order to calculate the masses of the scalar (a) and pseudoscalar (tt) bosons, we include the pseudoscalar degree of 
freedom in V and take the second derivative. The values of masses are obtained by multiplying the second derivative 
by the appropriate factor /: 



M 2 = 

a 3 da 2 



(56) 



Here, "min" at the end of the equations means that they are evaluated at the minimum of V(a, 05; mgl In this paper, 
we do not examine the factor /; rather, we study the "mass ratio" M a /M^. As discussed in Refs. @Q, below T c , 
M,r(T) weakly depends on the temperature and the value is dominated by the current quark mass and for T > T c 
the pion loses its Goldstone nature. On the other hand, M a (T) decreases in association with the partial restoration 
of chiral symmetry. Therefore, above some temperature T* , M a (T) becomes smaller than 2M V (T) and the width 
r<j^27r \J\ — AM 2 /M 2 vanishes. Figure 9 shows the temperature dependence of M a /M n . We evaluate T* and find 

T* ~ 0.97T C . (57) 

If we fix the Af 7r (T = 0) to 140 MeV, then, M a (T = 0) turns out to be 668 MeV. The temperature dependence of M a 
and M„ is shown in Fig. 10. 



B. T = 0,fi ^ case 

Figure 11 shows the chemical potential dependence of the effective potential at T — for mfl(lGeV) = 7 MeV. 
In the chiral limit, the critical chemical potential is (i c = 422 MeV. We note that, in m,R(lGeV) = 7 MeV case, 
another extremum appears for relatively larger value of (i. Figure 12 shows the (i dependence of cr min . The quark 
condensate discontinuously vanishes at (i — (i c in the chiral limit, indicating the first order phase transition. In 
mfl(lGeV) = 7 MeV case, a m i n changes discontinuously at (i = /i*, a little larger value than (i c . However, because of 
finite quark mass, it does not vanish for (i > (i*. It is widely expected that as rriR, namely the external field, grows 
the discontinuity weakens or disappears. We examined up to ttir ~ 400 MeV, however, the discontinuity does not 
vanish. 



C. r/0,/i/0 case 

Figure 13 shows the (i dependence of tr TO j„ for several temperatures. We expect that there is the critical end point 
E, where the discontinuous jump of the order parameter ends in the T-(i plane. Apparently its position is 

T E - 100 MeV, (i E ~ 300 MeV, (58) 

though it is not easy to determine accurately. The physics near the critical end point is interesting in the light of 
relativistic heavy-ion collision experiments pj. This, however, is beyond the scope of this work. 
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IV. SUMMARY AND DISCUSSION 



The primary purpose of this paper was to study the current quark mass effects on chiral phase transition of QCD 
in the improved ladder approximation. We made use of the CJT effective potential because the use of the SDE only 
is not adequate for studying the phase transition, and moreover, it was known that, in studies of the SDE, there is a 
difficulty in removing the perturbative contribution from the quark condensate |K^]i"7f ] . 

To begin with, we modified the form of the CJT effective potential using the two representations of the SDE. 
We then studied the tn and ttir dependence of the order parameter and obtained the reasonable results. Being 
motivated by Ref. [^J, the critical number of massless flavors was also studied. Our result is consistent with lattice 
QCD results for relatively small £r. Incidentally, our formulation of the effective potential is entirely based on the 
Higashijima-Miransky approximation. However it was known that the chiral Ward-Takahashi identity is broken unless 
gluon momentum squared is used as the argument of the running coupling Therefore, it is desirable to formulate 
the effective potential for the finite current quark mass independently of the Higashijima-Miransky approximation. 

We then extended the effective potential to finite temperature and density by introducing the functions 2?t,/x(p), 
£t>(p) and S 5; t, m (p)- We calculated the effective potential numerically and investigated current quark mass effects 
on the chiral phase structure. In either T^O, /i = 0orT = 0, /j^O case, the behavior of the order parameter in 
low-T or \ow-fj, region does not suffer the effects of finite, however small enough, current quark mass. On the other 
hand, in high-T or high-/x region its behavior is perfectly dominated by the current quark mass. We also examined 
numerically the critical exponents /3, 7 and 6 and confirmed the Landau theory of the second-order phase transition. 
The temperature dependence of the mass ratio M a jM^ was also studied. We found that for T > T* = 0.97T C the 
mass ratio becomes smaller than two. In Ref. [ p5|]38| , the temperature T* has been obtained and our results are 
almost coincident with theirs. We note that the similar temperature dependence of the mass ratio have been also 
obtained from the Nambu-Jona-Lasinio model p4| ] and the linear sigma model p^ ] . In the previous paper |l5| ] , in the 
chiral limit, we found the tricritical point at T P = 107 MeV, ^ = 210 MeV. We found ,in m fl (lGeV) = 7 MeV case, 
the position of the critical end point is Te ~ 100 MeV, /tig ~ 300 MeV. 

Finally, some comments are in order. The treatment of the quark and gluon propagator at finite temperature and 
density is somewhat oversimplified in the present work. In the future, we would like to consider the wave function 
renormalization and more appropriate functional form for Y,t,^(p), X5 : t,^(p) and T>t^{p) which would depend on 
T and fi explicitly. In particular, we should take into account the screening of the gluon in T>t^{p), that is to say, 
the effects of the Debye (electric) mass and the nonperturbative magnetic mass, which arises from the dimensional 
reduction, m mag ~ g 2 T at T ^ [p6| . They may affect the precise location of the critical line. We also note that 
these screening masses are probably related to the infrared regularization parameter Pr, though it is not easy to show 
it explicitly. Furthermore, it is well known that gauge covariance is lost, namely the Ward-Takahashi identity for the 
fermion-gauge-boson vertex is not satisfied when one uses the ladder approximation. Of course, gauge invariance in 
the case of the non-Abelian gauge theory is ensured by the Slavnov- Taylor identity which, in covariant gauges, includes 
a contribution connected with the Faddeev-Popov ghost fields. Thus, addressing ourselves to gauge covariance in QCD 
is difficult, however, progress has been made with the Abelian gauge theory jL6|. Although we suppose that results of 
the analysis will be essentially unchanged, it is preferable to use the correct form of the fermion-gauge-boson vertex. 
As concerns the phase structure, it is interesting to study the physics near the critical end point, for example, critical 
slowing down J37J. We also plan to study the quark pairing including a color superconductivity and a "color-flavor 
locking" § (for N c = N f = 3 case). 
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APPENDIX: 

In this appendix, we show the effective potential V(a, 05; uir) explicitly. In the first place, we consider the case of 
zero temperature and finite chemical potential. 
Using Eqs. (|^) and |l|), we obtain 
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Vi(a, a 5 ;m R ) 



d 4 PE St, p (p) + K.t,M + (P4 + *M) 2 + lif 
(2tt)* 111 



1 

47T 3 

-<5« 



In 



(p 4 + i[i) 2 + |p| 2 
(g^fr) + E 2 :T » + pj + |p1 2 - ^ 2 ) 2 + (2 m ) 2 

(p2 + bT - m 2 ) 2 + (2m) 2 



(Al) 



(A2) 



Moreover, in order to remove the divergence, we introduce 5^ = Vx(0, 0;?tir) the term independent of a and cr.5. 
In Eq. (|37|), we carry out the momentum differentiation and, then, use the Eqs. (fi3|), ( ^i|) and (|45|). V 2 is obtained 

as 



where the imaginary part of V\ is odd function of p^\ therefore it has been removed from Eq. ( |Al[ ) and 

/•oo 

dp 4 / dip] bf . 

) Jo 



V 2 (a, <T 5 ;m R ) 



(pi + \p\ 2 ) 2 (pi + bT +P 2 7 )[ln(p! + 2 

2 i U2 x „2 1 i^„2 j_ |„12 i „2 \ i J i l„12 ° \P' a i m R) 



3C 2 air s J p (pj + \p\ 2 + p 2 r ) Hp 2 + \p\ 2 + p 2 R ) +pl + \p\ 
(p 2 + b1 2 ) 2 [ln(p 2 + |p1 2 +^)r 2 



4-1 



3Tr 3 C 2 a J p {p 2 + \p\ 2 + p 2 R ) Hp 2 + \p\ 2 + p%) +pl + W 
[Hp 2 + \p\ 2 +p 2 R ) + l-a/2] 2 



(pI + \p\ 2 +p r ) 3 



(A3) 



where the function 

S{p;a,m R ) 



am r 



[HpI + \p\ 2 +p 2 r )} 



-a/2-1 



pI + \p\ 2 +p r 



{pI + W+pr) 2 



[Hpl + \pf+p 2 R )\ 



2 00/2-2 



ln(p 2 + |pl 2 +p 2 J ) + l 



(A4) 



is introduced and (5^ 2 ) = V 2 {Q 1 0; tur) is again the subtraction term indepen dent of a a nd a 5. 

At finite temperature and chemical potential, the p 4 integration in Eqs. ( |AlD and (A3) is replaced by a sum over 
the fermion Matsubara frequency. 
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FIG. 1. The t R dependence of — (qq) 1 ^ 3 in the chiral limit. The critical value of t R is (t R ) cri ucai = 0.56. Notice that — {qq} 1 ^ 3 
is taken to be dimensionless. 
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FIG. 2. The tur dependence of (qq) for the cases t R = 0.1 (solid line) and t R = 0.8 (dotted line). All the quantities are 
taken to be dimensionless. 
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FIG. 3. The in dependence of n°f. shows a tendency to decrease as tn grows. 
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FIG. 4. Schematic view of the phase diagram obtained from our model in the chiral limit [15]. We have found a second-order 
transition at T c — 129 MeV for /i = and a first-order one at fi c — 422 MeV for T = 0. Solid line indicates the phase transition 
of second-order and dashed line indicates that of first-order. A tricritical point P has been found at Tp — 107 MeV, pip = 210 
MeV. 
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FIG. 5. The effective potential at finite temperature and zero chemical potential. V is defined by V = 24-7r 3 V. The curves 
show the cases T/A QCD = 0, 0.15 and 0.2 with m R {lGeV) = 7 MeV. 




FIG. 6. The temperature dependence of a m in at /x = 0. The curves show the cases mfl(lGeV) = 7 MeV (solid line) and the 
chiral limit (dotted line). 
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FIG. 7. Linear log fit to ln(qq) \m R =o (circle) and lna2(T) (plus). Here the absolute normalization of (qq)\m R =o and flS2(T) are 
taken arbitrary. The gradient of the linear function fitted the data for ln(gg) | mH =o and lnci2(T) are /3 = 0.4997 and 7 = 1.000, 
respectively. 
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FIG. 8. Linear log fit to ln(qq) |t=t c • The absolute normalization is taken arbitrary. The gradient of the linear function 
fitted the data is 1/6 = 0.3329. 
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FIG. 9. The temperature dependence of the ratio M a /M n - 




FIG. 10. The temperature dependence of M CT (solid line) and M n (dotted line). 
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FIG. 11. The effective potential at finite chemical potential and zero temperature, where V = 2Att 3 V. The curves show the 
cases ;u/Aqcd = 0, 0.5 and 0.6 with m R (lGeV) = 7 MeV. 



0.5 




J 0.2 - 



0.1 - 

' 1 1 1 1 — 

100 200 300 400 500 

H [MeV] 

FIG. 12. The chemical potential dependence of a min at T = 0. The curves show the cases mfl(lGeV) = 7 MeV (solid line) 
and the chiral limit (dotted line) . 
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